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Abstract 

The mass power spectrum for a Universe dominated by the Chaplygin gas is 
evaluated numerically from scales of the order of the Hubble horizon to 100 Mpc. 
The results are compared with a pure baryonic Universe and a cosmological constant 
model. In all three cases, the spectrum increases with k, the wavenumber of the 
perturbations. The slope of the spectrum is higher for the baryonic model and 
smaller for the cosmological constant model, the Chaplygin gas interpolating these 
two models. The results are analyzed in terms of the sound velocity of the Chaplygin 
gas and the moment the Universe begins to accelerate. 

-xf : 

^ ■ PACS number(s): 98.80.Bp, 98.65.Dx 

o : 

§ ■ 1 Introduction 

Since the observations of the supernova type la have indicate that the Universe must be in 
an accelerated expansion today [l], |2| the nature of the fluid responsible for this inflationary 
behaviour has been object of many studies. In order this accelerated expansion to take 
place today, the Universe must be dominated by a fluid of negative pressure. Moreover, 
.£h ! this fluid must remain a smooth component of the matter content of the Universe, since 
it does not appear in the dynamics of clusters of galaxies, for example. The most natural 
candidate for this "dark energy" component is a cosmological constant [Q]. However, there 
are two main problems with the idea that a cosmological term must be the dominant 
component of the Universe today: First, in spite of the fact that quantum field theory 
predicts a cosmological constant, the observational value would be 120 orders of magnitude 
smaller than the predicted value; second, it is a quite special fact that the energy density 
associated to a cosmological constant have a value near to the density of other types of 
matter exactly today, a fact that is usually called " cosmic coincidence" . A very popular 
way of coping with these problems is through the introduction of a self-interacting scalar 
field, called "quintessence" J|, ||. For some kind of potential terms, which have their 
justification in fundamental theories like supergravity, there are the so-called tracking 
solutions: The self-interacting scalar field evolves in such a way that it approaches a 
cosmological constant behaviour exactly today ||. But the great variety of possible 
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potential terms with these properties, and the necessity of fine tuning some parameters 
in order to generate them, makes quintessence a promising scenario, but without any 
definitive model until now. 

In a previous work J7|, we have considered the possibility that the Universe today is 
dominated by another kind of fluid, called Chaplygin gas. The Chaplygin gas is charac- 
terized by an equation of state 

p = , A = constant . (1) 

P 

Hence, it exhibits negative pressure, as required to explain the acceleration of the Universe 
today. However, the fact that it depends on the inverse of the density gives it some special 
properties. In particular, the conservation of the energy-momentum tensor implies that 
the Chaplygin gas density depends on the scale factor of the Universe a as 




P=\ A +~« > (2) 



where B is an integration constant. This relation, first presented in M where some cos- 
mological consequences were explored, has some interesting properties: For small values 
of the scale factor, the Chaplygin gas exhibits the same behaviour as a pressurelless fluid; 
for large value of the scale factor, it approaches the equation of state of a cosmological 
constant. Hence, in some sense it reproduces the quintessence scenario, interpolating a 
matter dominated Universe and a cosmological constant dominated Universe, but in quite 
different way, since no appeal is in principle made for a self-interacting scalar field. 

The Chaplygin gas has been first identified in the study of adiabatic fluids B. More 



recently, it has received an interesting motivation, connected with string theories |]10|, |Tl 
Considering a D-brane in a D + 2-dimensional space-time and employing the light-cone 
parametrization, the Nambu-Goto action reduces itself to the action of a newtonian fluid 
which obeys the equation of state (|I|). Even if this reduced action is that of a newtonian 
fluid, the related group of symmetry has the same dimension as the Poincare group, 



revealing that the relativistic origin is somehow hidden in that equation of state [|12||. In 



this sense, when we speak about the Chaplygin gas, we can be speaking about a gas of 
.D-branes inaD + 2-dimensional space-time. 

The interest of this equation of state for cosmology can be summarized by evoking two 
main properties of the Chaplygin gas in an expanding Universe. First, as it has already 
been said, it behaves initially as pressurelles fluid, and later as a cosmological constant. In 
this sense, it plays the same role as the scalar field in the quintessence program. Second, 
the sound velocity of the Chaplygin gas is positive, in spite of the fact that its pressure is 



negative. This a very important property, since, as it has already been shown in [13|, fluids 
with negative pressure obeying a barotropic equation of state suffer from instabilities at 
small scales due to an imaginary sound velocity. 

In [0, the evolution of density perturbation in a Universe dominated by the Chaplygin 
gas has been analyzed exploiting the fact that the Chaplygin gas obeys the newtonian 
equations of hydrodynamic. It has been shown that the predicted scenario is compatible 
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with the formation of local structures, and that the perturbations in the Chaplygin gas 
tends to vanish asymptotically, as we should expect due to the asymptotical behaviour 
similar to a cosmological constant. In the present paper, we intend to go a little further. 
Employing the relativistic equation for the description of the evolution of the Universe, 
considering a matter content composed of cold dark matter, baryons, radiation and the 
Chaplygin gas, we study the evolution of the power spectrum of the clustered matter. 
Starting with a given primordial power spectrum, which is very close to the Harrison- 
Zeldovich spectrum, at the end of the radiative phase, we show that this power spectrum 
evolves to an increasing one at moments near today. In order to track the specificity of role 
played by the Chaplygin gas, we perform the same analysis to the case where the Chaply- 
gin gas is replaced by a cosmological constant and to the case of a Universe composed of 
radiation plus pressurelles matter only. The results indicate that the spectrum predicted 
by the Chaplygin gas model interpolates between that of a baryonic Universe (which pre- 
dicts a strong increasing spectrum today for the scales analyzed here and the primordial 
spectrum employed) and a cosmological constant model (which predicts a moderately 
spectrum today for the same scales) as the sound velocity associated with the Chaplygin 
gas goes from zero (A = 0) to the velocity of light. 

The model developed here has some important simplifications. We consider, for ex- 
ample, the neutrino and photon contribution to the radiative fluid, whithout taking into 
account the different moments of decoupling of these different components. Moreover, we 
do not take into account that after the beginning of the material phase, photons remain 
coupled to the baryonic matter until the decoupling moment. In the same way we do 
not take into account that the cold dark matter fluid decouples from radiation before the 
baryonic fluid. We restrict our analysis to scales that goes from 100 Mpc until the Hubble 
horizon at the present moment, and we do not consider the neutrino free streaming and 
photon diffusion process that occurs when the perturbations are inside the horizon. All 
these simplications limit the possibility of comparing in detail the obtained results with 
the recent surveys of clustered matter. But we emphasize that the goal of the present 
paper is to identify some main features of the Chaplygin gas concerning the structure 
formation problem, and in view of this goal those simplifications (even if very strong from 
the observational point of view) do not spoil our main conclusions. However, the results 
are qualitative in agreement with more precise and realistic calculations. 

This paper is organized as follows. In next section, we describe the cosmological 
model with the Chaplygin gas, comparing it with a baryonic Universe and a cosmological 
constant model. In section 3, we settle out the perturbative equations which are analyzed 
numerically. The power spectrum is evaluated for each of the models described before. 
In section 4, we present our conclusions. 

2 The background model and the cosmological pa- 
rameters 

Let us consider a quite general model consisting of a pressurelles fluid, radiation, a cosmo- 
logical constant term and the Chaplygin gas. If we take a Friedmann-Robertson- Walker 
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metric, which describes an homogeneous and isotropic Universe, the Einstein's equations 
reduce to 



a\ 2 k 8nG , . , . 

- +— = —z-{Pm + Pr + PA + Pc) , (3) 
ay cr 6 

- + 2(-) +\ = -87rG(p m +p r - p A +p c ) , (4) 

where a is the scale factor of the Universe, k is the curvature of the spatial section, and 
p m , p r and p c are the density for the pressurelles fluid, radiation and the Chaplygin gas, 
p m = 0, p r = p r /3 and p c = —A/p c being their corresponding pressures, pa is the density 
associated to the cosmological term. In what follows we will be interested in three cases: 
Pa = 0, the Chaplygin gas model; p c = 0, the cosmological constant model; p\ = p c = 0, 
the baryonic model. 

The curvature of the spatial section will be set equal zero. This simplification is 
justified by the fact that the data concerning the anisotropy of the cosmic microwave 
background radiation favor a flat, or almost flat, Universe. In fact, the parameter associ- 
ated with the curvature of the spatial section of the Universe is defined as f2fc = 1 — Qt, 
where Qt is the total observed density of the Universe divided by the critical density. 
The position of the first Doppler peak of the spectrum of the anisotropy of the cosmic 
microwave background radiation is directly connected to Qt- The recent results from the 
BOOMERANG and MAXIMA measurements of CMB spectrum give tt k = ± 0.06 fT?] . 



Hence, we can fix k = without to oversimplify our model. Moreover, an inflationary 
phase in the primordial Universe predicts (except for very special cases) k = 0. 

The different fluids are assumed to interact through the geometry only. Hence, the 
energy-momentum tensor for each component conserve separately, 



7T ;M = , (5) 



leading to the expressions 




p m = ^rr ' Pr = ^T i Pc = \IA+— , p A = constant . (6) 
a 6 a 4 V a 

Henceforth, we will fix the value of the scale factor today as equal to one, ao = 1. In this 
case, pmo, Pro and p c o = \J A + B are the density for pressurelles matter, radiation and the 
Chaplygin gas today. In the expression for the Chaplygin gas, the integration constant B 
can be expressed in terms of A and p c0 , leading to B — p 2 c0 — A, what allows to rewritte 
the expression for the Chaplygin gas density as 

p c = p c0 \IA + ^^ , A = 4r ■ ( 7 ) 

PcO 

The redefined parameter A will be connected with the sound velocity for the Chaplygin 
gas today as we will see later. 
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Using the equations of motion and the previous result, the first and the second deriva- 
tives of the scale factor can be expressed in terms of the mass parameters and the scale 
factor itself: 



F(a) =<R 



Q-mO &<r0 ,2— - ' ' ~ ^ ~ A 



^ + a 2 n A0 + n c0 \ a 4 A 



a a 2 v a z 



1/2 



G(a) = ^{-^-2^f + Q c0 2M3 ^ +2n A0 a\ , (9) 
- ' a « 3 - Ma 4 + ±=£ 



where Q r Q, Q m o, f2 c o and Q A q are fractions of the radiation density, pressurelles matter 
density, Chaplygin gas density and cosmological constant density with respect to the total 
mass today, respect ivelly. Of course, fl r0 + Q m0 + Q c0 + Qa = 1, since we have assumed 
a flat Universe. Moreover, R = 87rGtgp /3, where t is the present age of the Universe 
and p is the total density of the Universe. In this way, (^],|9]) have become dimensionless 
equations, the dots meaning derivatives with respect to t/t , that is, fractions of the age 
of the Universe or, if one prefers, the redshift z. 

In general, expressions (f||9|) do not admit simple closed solutions, except for the case 
where Q c q = Q A q = 0. However, since in more general cases numerical integration will be 
needed, we will not consider that particular case, which, in any way, leads to background 
solutions that do not allow simple analytical integration of the perturbed equations. How- 
ever, expressions will be crucial in settling out the perturbed equations which will 
be used in the numerical integration procedure. 

Before to treat perturbatively the system described before, the cosmological param- 
eters for each model will be specified. First of all, the Hubble parameter today will be 
taken as 

km , . 

#o = 72- . (10) 

Mpc.s 

This leads to an age of the Universe of the order 

t ~ 4.1 x 10 17 s ~ 13.3 Gy . (11) 

This value is quite favored by the recent estimations of the anisotropy of cosmic microwave 
background, supernova type la observations, being also in good agreement with recent 
re-estimations of the age of galactic objects obtained through the Hipparcos program 



Tq, 19 1 . Using the Einstein's equation, the resulting total density of the Universe is 



Po ~ 10 -29 g/cm 3 . In the next section, we will analyze three main models: 

• Model I: Q c q = Q A q = 0, which will be called baryonic model; 

• Model II: ^ao = 0, which will be called Chaplygin gas model; 

• Model III: Q c q = 0, which will be called cosmological constant or ACDM model. 

Using estimations of cosmological parameters, coming from the various observational pro- 
grams in course, as summarized in [|T7[], we will fix the following values for each model: 
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• Model I: fi r0 = 10" 4 ; fi m0 = 0.9999; 

• Model II: Q r0 = 10~ 4 ; Q m0 = 0.2999; Q c0 = 0.7; 

• Model III: Q r o = 10" 4 ; n m0 = 0.2999; Q A = 0.7. 

The values used above are not the unique possibility of course. But, they are consistent 
with the different observational estimations: The data used above are well inside the 
observed values, taking into account the bar errors. For the estimation of the radiative 
contribution, neutrino and photon components are taken into account. For the pressurel- 
less fluid, baryonic and weakly interacting massive particles (WIMPS), which must be 
the components of the cold dark matter, are considered. In model II, specifically, the 
possible contribution of the cosmological constant is replaced by the Chaplygin gas, while 
in model I both exotic fluids are ignored, only radiation and pressurelless matter being 
kept. 

Each of the models described before imply a different moment for the equilibrium 
between radiation and the other forms of matter. For a pure cosmological constant model, 
this equilibrium occurs for z eq ~ 3, 000. For the Chaplygin gas models, this moment 
depends on the value of A: For A = 0.5, it occurs at z eq ~ 4, 600, for A = 0.8 at 
z eq ~ 6, 100, and for A = 0.95 at z eq ~ 8, 000. The equilibrium between radiation and 
matter is reached at z eq ~ 10, 000 for the baryonic model. These different values for z eq 
do not play a significant role since the background models imply a smooth transition from 
the radiation dominated era to the matter dominated era, the transition extending for 
a non negligible period of time. Hence, the exact moment where the radiation density 
equals the matter density is of minor importance in this scenario. 

An important point must be noticed. In model I, there is an initial radiative phase, 
followed by a material phase. The Universe remains always in a subluminal expansion. 
On the other hand, for models II and III, there is a moment, quite near today, where 
the subluminal expansion gives place to a superluminal expansion, that is, inflation. This 
moment can be easily determined putting G(a) = 0. For the cosmological constant model, 
with the values given above, inflation begins to occur at z ~ 0.67. For the Chaplygin 
gas model, this depends crucially on the value of A: For A = 0.5, it occurs at z ~ 0.01; 
for A = 0.8, it occurs at z ~ 0.22; for A = 0.95, it occurs at z ~ 0.43. The greater A, 
the earlier the superluminal expansion of the Universe begins. This is natural, since as 
A — > 0, the Chaplygin model coincides with the baryonic model, where no superluminal 
expansion occurs at any time; on the other hand, as A — ► 1, the Chaplygin gas model 
coincides with the cosmological constant model. 

Finally, in what concerns specifically the Chaplygin gas model, there is an initial 
radiative phase, followed by a matter dominated phase. The final stage is a deSitter 
phase, with the Chaplygin gas playing the role of a cosmological constant of value p\ = 
\J~A~. The Chaplygin gas behaves initially as a dust fluid, but it implies an intermediate 
phase characterized by a mixing of stiff matter and cosmological constant [Ef. In ET 



a "generalized Chaplygin gas" was studied, considering an equation of state of the type 
p = —A/(p a ), where that intermediate stage becomes a mix of cosmological constant and 
a fluid characterized by p = a p. In spite of being interesting to consider this more general 



6 



equation of state, we will restricted ourselves to the traditional Chaplygin gas, with an 
equation of state given by (|l|). 



3 Perturbative analysis and the mass power spec- 
trum 

In order to perform the perturbative analysis necessary to determine the mass power 
spectrum today, we return back to the einsteinian equations, writing them as 

R, l/ = 8nGj2(T; u -lg, u T) , T^ = , (12) 

where the indice i labels the different matter components described in the preceding 
section. The background solutions are described by the functions The perturbations 

are introduced by writting the dynamical variables as 

o o o 

SV = 9nu +h llu , Pi = Pi +5 Pi , Pi = Pi +Spi , (13) 



where g^, Pi and Pi represent the background solutions and h^ u , dpi and dpi are small 
fluctuations around them. Due to the coordinate reparametrization freedom, we can 
impose a coordinate condition, which we choose to be the synchronous one: h^o = 0. 

We will consider just adiabatic perturbations. In this case, the velocity of the sound 
for each component is given by 

v h% • <") 

For the Chaplygin gas, this reads 

4-4 . as) 

& pi 

which is positive. Hence, the sound velocity of the Chaplygin gas is given today by 
v s = vAc. As A approaches one, the sound velocity of the Chaplygin gas approaches 
the velocity of light and, by (|7|), the Chaplygin gas becomes essentially a cosmological 
constant. 

After a long but standard calculation JTJ], [15|, [HJ, we end up with the following system 



of coupled perturbed equations for the general model composed of radiation, pressurelless 
fluid, Chaplygin gas and cosmological constant: 



h + = 3Q 2 j2A r a. + A m n m +(l + 3^A c fi c | , (16) 



A, + ^(e r -|) = , (17) 

e r + ^e r = , (is) 

a 4 cr 
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h 

2 



A 



e,-- 



2 + 34)^e e 



a A 
-6--A c 

a P 2 C 
A 



p 2 c - A a 2 



(19) 
(20) 
(21) 



In these equations h = Aj = — is the density contrast for each fluid component, Q{ = 
denotes the mass fraction for each fluid at a given time, being the total mass density 
at that time, Bj is the perturbed velocity for each fluid. A Fourier decomposition has also 
been performed, writting all perturbed quantities as Sf(x,t) = 8f(t)e~ ik ' x , k being the 
wavevector associated with each Fourier mode. Notice that the fluid velocity perturbation 
for the pressurelless fluid and the density perturbation of the cosmological constant are 
zero. In particular, the cosmological constant influences the perturbed equations through 
the behaviour of the background quantities only. 

Since we do not have in the general case an analytical expression for the scale factor 
as function of the time, we rewritte the above equations in terms of the scale factor 
itself, using (§,0). Moreover, we make equations (16,17,18,19,20,21) dimensionless by 
multiplying them conveniently by to, the age of the Universe, and reinserting, when is the 
case, c. Finally, we eliminate h, expressing it in terms of A m , using ([19]). The perturbed 
equations read then 



A" - 

in 



a: 



(G(a) 
3lF(a) 



-\x. 



— \2A r tt r + A 



'l + 3E(a))A c tt ( 



- A' 



K+ 



1 - E(a) 



F(a) 
2 + 3E(a 



a 



A' 



80ttA; 



-6- 



k 

E(a 



A r 



4a 2 F(a) 



■A, 



a 

80nk 
k 



E(a) 



a 2 F(a)(l -E(a)) 



where the primes mean derivative with respect to a and E(a) = Aa 6 /[Aa 6 + 
mass parameter are written as 



n 



rO 



F 2 (a)a? 



i/iO 



F 2 (a)a 



n c = n c0 \ a + 



1 - A 



F 2 (a) 



,(22) 
(23) 
(24) 
(25) 
(26) 
A)]. The 

(27) 



The wavenumber has been conveniently reparametrized in terms of a reference scale k = 
2ir/\ , where A = 100 Mpc. Notice that, since we have imposed the value of the scale 
factor today equal to one, the co-moving wavelengths are equal to the physical wavelengths 
today. 

Equations (^,^,|^,^,^) represent the complete equations for the Chaplygin gas 
model, with ^ao = in (0,0). The cosmological constant and the baryonic models are 
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obtained by fixing A c = O c = in the perturbed equations above. For the cosmological 
constant model, it must be fixed Q c0 = in the unperturbed equations, while for the 
baryonic model it must also be fixed Q^o — 0. 

The perturbed equations written above are completely consistent mathematically with 
the background model described by (|8|||). However, from the physical point of view, 
there are many simplifications which must be pointed out. First, what is called here 
"baryonic" matter is a general word for "pressurelles fluid", and it includes both the 
baryonic matter properly speaking and the cold dark component (WIMPS). Even if the 
cold dark matter component can also be represented by a zero pressure fluid at the stage of 
galaxy formation, it decouples from the photonic fluid earlier than the baryonic matter due 
to its very small cross section. This fact has consequences for the mass power spectrum; 
but we neglect it here, since it will affect the final results in the same way for the three 
models to be studied. 

The second important approximation concerns the neutrino free streaming. Of course, 
the density of neutrinos will affect the evolution of the background. However, due to its 
decoupling at quite early times, the neutrinos will not contribute to the radiation pressure 
as the photonic component does in the gravitational collapse process. Hence, in principle, 
this should be taken into account in the perturbative equations. Again, this effect will 
be neglected: The radiation density in the perturbed equations will be taken as in the 
background equations. As in the preceding case, since this effect will be neglected in the 
same way in all three models to be studied, this approximation does not spoil significantly 
the results from the point of view of identifying the main features of the Chaplygin gas 
compared with the pure baryonic Universe and, mainly, with the cosmological constant. 
Another simplification concerns the photon diffusion, which will not be considered also, 
since it leads to the breakdown of the perfect fluid approximation. 

Having established the main features of the perturbative analysis to be carried out, 
we turn to some important aspects of the problem. Equations ( P^| , p5| , pi| , Po] , pBD are com- 



plicated enough to prevent us to try to obtain analytical solutions, even in the long 
wavelength limit. Hence, we will integrate them numerically. This integration is carried 
out from z ~ 10 4 (approximatelly the moment of equilibrium between matter and radia- 
tion) until z = (today). It remains the important question of the initial conditions. The 
initial velocity perburbation Oj will be set equal to zero at that moment, as well as the 
derivative of the perturbations of the matter components, = 0. For the matter pertur- 
bations themselves, the most natural choice, and the more consistent with respect to the 
results of the anisotropy of the cosmic microwave background, is the Harrison- Zeldovich 
spectrum, which states that all perturbations have the same amplitude when they cross 
the horizon. The Harrison- Zeldovich spectrum implies that |15| 



{bk)H oc k~ 3/2 . (28) 

However, it comes out to be more convenient to impose the initial conditions at a constant 
time tj, instead to impose them at the horizon crossing moment, since this crossing occurs 
at different moments for different wavelengths. In order to obtain a spectrum near the 
Harrison-Zeldovich spectrum at the horizon crossing, the spectrum at a constant time 
defined before this crossing must be different from fl2"8"|). In [I5[] two phases, the radiative 
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and material ones, were considered and it came out that the spectrum at constant time 
should be written as 

(h)t^k l l 2 , (29) 

in order to obtain a Harrison- Zeldovich spectrum at horizon crossing. In our model, 
the situation is more complex since there are not two so well definite phases. However, 
( p5D remains a reasonable approximation, as it has been verified numerically. It must be 
stressed, in any case, that the goal of the computation performed in this work is not to 
test the primordial spectrum, but to identify the specific features of the Chaplygin gas. 
In figures 1, 2 and 3 the behaviour of the density contrast for matter, Chaplygin gas and 
radiation are shown for the Chaplygin gas model with A = 0.8 and k/ko = 0.5. 

We have evaluated numerically the power spectrum today for the dust component, 
represented by the baryonic and cold dark matter component, given their power spectrum 
at t — ti (equivalent to z = 10 4 ). We restrict ourselves to these components since they are 
essentially what is observed in clusters and super-clusters of galaxies. In figure 4 we display 
the results for P(k) = k 3 ^ 2 5k as function of k. The numerical integration was performed 
from scales of the order of the Hubble horizon today, X max ~ ~ 4,000 Mpc for the 
value of the Hubble parameter chosen, down to A m j n = 100 Mpc. The normalization of 
each spectrum is made by its value at \ max - Hence, each curve in figure 4 has a slightly 
different parametrization. However, what is more relevant is the slope of each curve. In 
this aspect we notice that all three models reveal an increasing spectrum. The slope of 
the spectrum is higher for the baryon model, becoming smaller for the Chaplygin gas 
model, approaching the cosmological constant model as A — > 1. Hence, the Chaplygin 
gas interpolates a pure baryonic model and a ACDM model as the velocity of the sound 
associated with it goes from zero to the velocity of light. In figure 5 we plot the same 
spectrum for In P 2 (k) = \n(k 3 5l) as function of In A. Notice, from this figure, that the 
spectrum for the cosmological constant tends to be more flat, what is compatible with the 
CMB measurements [fTTH , deviating significantly from this "almost Harrison-Zeldovich" 
spectrum as we approach the pure baryonic spectrum, making the velocity of sound of 
the Chaplygin gas going from v 2 = c 2 to v 2 = 0. 

The interpretation of these results makes appeal for some physical aspects of the 
problem. Notice that the three models contain the same amount of radiative fluid. As 
the perturbation enters in the horizon, the radiative fluid contributes to the pressure and 
may distort the initial spectrum. In principle, this occurs in the same way for the three 
models. A quite significant effect may come from the moment where the perturbation 
enters in the horizon, which is evidently different for each model. The condition for the 
horizon crossing is given by the relation k = aH/c. Using the parametrizations described 
before, this relation takes the form k/k Q = F(a)/(40). A perturbation with a scale of 
the order of 100 Mpc today has entered in the horizon at z ~ 1,402 for the baryonic 
model, while for the cosmological constant model it has entered at z ~ 2,771. For a 
perturbation of the order of 2, 000 Mpc today, the horizon crossing occured at z ~ 12 
for the cosmological constant model, and at z ~ 3 for the pure baryonic model. The 
Chaplygin gas models leads to values between those two. As the cosmological constant 
model is approached, the perturbations analyzed have entered earlier in the horizon, being 
subjected for a longer time to the effects of pressure. 
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Another important effect, which seems to be more significant than the previous one, 
comes from the moment the Universe begins to accelerate. The expansion of the Universe 
leads to a friction like term in the perturbative equations (see, for example, (|T6|)). The 
accelerated expansion tends to smooth the fluctutations more efficiently than a non- 
accelerated expansion. The baryonic model never accelerates, as it was discussed in the 
preceding section. The cosmological constant model begins to accelerate quite recently, 
but in such a way that the fluctuations in density are significantly affected. Of course, the 
Chaplygin gas models leads to an accelerated Universe whose time transition is near to 
our present time if A is small, and near to the cosmological constant model time transtion 
as A approaches one. This fact is reflected in the final evaluated spectrum, in special in 
the slope of the spectrum. 

4 Conclusions 

The Chaplygin gas appears as an alternative to the cosmological constant and to quintessence 
in order to explain the stage of accelerated expansion of the Universe today. It is a fluid 
with negative pressure which depends on the inverse of the density. Due to its equation 
of state, it interpolates a matter dominated phase and a cosmological constant phase in 
what concerns the evolution of the Universe. In this paper we have addressed the ques- 
tion of structure formation in the realm of a Universe dominated by the Chaplygin gas. 
The goal of the present work was to determine the features of the mass power spectrum 
inside the horizon today given a primordial spectrum in a given time near the equilibrium 
between radiation and matter, around z ~ 10 4 . In order to track the specific behaviour 
of the Chaplygin gas, we computed also the mass spectrum for a pure baryonic Universe 
and for a ACDM model. 

An important parameter in the Chaplygin gas model is played by the constant A in 
the equation of state of this fluid. It is related to the sound velocity in such a fluid. 
The sound velocity of the Chaplygin gas today is given by A = where p c o is the 

PcO 2 

Chaplygin gas density today. If A = 0, the Chaplygin gas behaves always as a dust fluid; 
as A approaches one, the velocity of the sound approaches the velocity of light, and the 
Chaplygin gas become indistinguishable from a pure cosmological constant. 

This fact is reflected in the mass power spectrum. We have evaluated this spectrum 
numerically, using an almost initial Harrison- Zeldovich spectrum. Considering scales be- 
tween the Hubble horizon and 100 Mpc (the normal threshold scale above which observa- 
tions indicate that the Universe may be considered essentially homogeneous and isotropic), 
we find that, for the models studied here, the spectrum increases with the wavenumber 
k. However, the slope is higher for a pure baryonic Universe, and smaller for a ACDM 
model. The Chaplygin gas model interpolates between those two as A goes from zero to 
one. This smoothing of the spectrum can be understood by remembering that the expan- 
sion of the Universe plays the role of a friction term, and an accelerated expansion tends 
to suppress more efficiently the local irregularities than a non-accelerated expansion. For 
the chosen values of the cosmological parameters, in special the mass parameters, the 
more one approaches a cosmological constant model, the earlier the Universe begins to 
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accelerate, and consequently the power spectrum becomes more flat. We notice, more- 
over, that the earlier the Universe begins to accelerate, the earlier the perturbation of a 
given scale enters in the horizon, being subjected no longer only to the the expansion but 
also to pressure. 

A next step would be to compare the results obtained for the spectrum with the obser- 
vational data. For the scales chosen, this means that high redshift surveys of structures 
must be compared with the theoretical models. Important programs in this sense are in 
course p0| , f2~!| |23|| . The data concerning the mass power spectrum are resctricted mainly 
to scales up to A ~ 800 Mpc [^], but the uncertainty is higher the higher is the redshift. 
This observational data cover only the higher values of the wavenumber in the spectrum 
we have determined here. But, we think that no detailed comparison, even in this in- 
terval, of our results with the observational data may be possible until a more realistic 
perturbative analyis is performed, mainly in what concerns to take into account the spe- 
cific features of the cold dark matter WIMPS and the neutrino free streaming. In very 
small scales, non-linear effects must also be taken into account. Notice that, in spite of 
all the simplifications, the qualitative features of the spectrum found are similar to the 



observational results (see ||22| , for example, and also [24]). In this sense, the present work 
is just a first step in the investigation if a Chaplygin gas model can be competitive with 
other models, like ACDM or quintessence. 
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Figures captions 



Figure 1: Behaviour of the density contrast for the ordinary matter, with A = 0.5 
and k/ko = 0.5. 

Figure 2: Behaviour of the density contrast for the Chaplygin gas, with A = 0.5 
and k/ko — 0.5. 

Figure 3: Behaviour of the density contrast for radiation, with A = 0.5 and k/ko = 



0.5. 



Figure 4 



The power spectrum P(k) = k 3 ^ 2 5k as function of k/ko. 



Figure 5 



Behaviour of In P 2 



(k) as function of lnA/Ao- 
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